J/ty production in ultraperipheral Pb+Pb and p+Pb collisions at LHC energies 
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We compute cross sections for incoherent and coherent diffractive J/vp production in ultrape- 
ripheral nucleus-nucleus and proton-nucleus collisions using two different dipole models fitted to 
HERA data. We obtain a reasonably good description of the available ALICE data for coherent 
J/ty production and present our prediction for the incoherent cross section. We also find that while 
the normalization of the cross section depends quite strongly on the dipole model and vector meson 
wave function used, the rapidity dependence is very well constrained. 
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I. INTRODUCTION 

The Color Glass Condensate (CGC) provides a con- 
venient way to describe strongly interacting systems in 
the high energy limit, where nonlinear phenomena, such 
as gluon recombination, become important. Because the 
gluon density scales as ~ A 1 / 3 , these nonlinearities are 
enhanced when the target is changed from a proton to a 
heavy nucleus. 

The structure of a hadron can be studied accurately in 
Deep Inelastic Scattering (DIS) where a (virtual) pho- 
ton scatters off the hadron. A large amount of pre- 
cise high energy electron-proton data measured at HERA 
has shown that the gluon density inside a proton grows 
rapidly at small Bjorken x or, equivalently, at high en- 
ergy. These accurate measurements have also been a cru- 
cial test for the CGC, and recent analyses have confirmed 
that the CGC description is consistent with all the avail- 
able small- a; DIS data [1, 2]. 

In order to apply the CGC formalism in the heavy 
ion environment one would ideally want to study DIS off 
heavy nuclei at high energy. The proposed LHeC [3] and 
eRHIC [4] experiments aim to perform these measure- 
ments. Before that, one can hope to obtain information 
about the dense gluonic matter in the nucleus by study- 
ing e.g. single [5, 6] and double inclusive [7-9] particle 
production in proton- nucleus collisions. These hadronic 
processes are, however, not ideal for precise studies as 
the parton level kinematics is not fixed by the final state 
particles, in contrast to DIS. 

Ultraperipheral heavy ion collisions, where two heavy 
nuclei barely touch each other, offer an interesting pos- 
sibility to study photon-nucleus scattering without a 
lepton-ion collider. Recently, the ALICE collaboration 
has measured [10] diffractive vector meson production 
in these collisions, opening a new possibility to study the 
small- a; structure of heavy nuclei. Diffractive events are 
especially interesting as they allow one to also study the 
transverse spatial distribution of the gluons in the nu- 
clear wave function, if the momentum transfer t can be 
measured. 

In this brief note we present our predictions for coher- 
ent and incoherent J/^ production in heavy ion and pA 



collisions. In Sec. II we first review the dipole models 
used here to compute photon-proton scattering. In Sec. 
Ill we recall the results for coherent and incoherent vec- 
tor meson production in 7A collisions and discuss how to 
compute the full nucleus-nucleus cross section. Finally in 
Sec. IV we present our results before concluding in Sec. 
V. 



II. DIPOLE CROSS SECTIONS 

In ultraperipheral nucleus-nucleus collisions strong in- 
teractions are heavily suppressed, and the electromag- 
netic interaction is expected to dominate. One can thus 
consider one of the nuclei as a source of (virtual) photons 
that scatter off the other nucleus. In the dipole picture 
photon-nucleus scattering is described as the photon fluc- 
tuating into a quark-antiquark color dipole which then 
scatters off the target nucleus. The dipole model is valid 
only when the Bjorken x of the gluon is small, we imple- 
ment this constraint by computing the cross section only 
when x < 0.02. 

In the literature there are many parametrizations avail- 
able for the dipole-proton cross section 

dcTj; 

"ToT" O^Ti i*Tj %p) — 2A/ (r^, bp, xp), (1) 

where Af is the imaginary part of the forward dipole- 
proton scattering amplitude, tt is the transverse size of 
the dipole, hx is the impact parameter of the 7-p collision 
and xp is the usual Bjorken variable of DIS in a diffrac- 
tive event. The dipole amplitude JV satisfies the BK [11] 
evolution equation, and ideally one would want to fit the 
initial condition of the BK evolution to the available DIS 
data (as done in Ref. [1]), solve the BK equation and use 
the obtained dipole amplitude when computing other ob- 
servables, such as diffractive vector meson production. 

However, computing diffractive events requires knowl- 
edge about the impact parameter dependence of the 
dipole amplitude. Straightforwardly including impact 
parameter dependence into the BK equation leads to an 
unphysical growth of the size of the proton with the evo- 
lution [12] unless this is regulated by hand at the confine- 
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ment scale [13, 14]. Because of this complication we use 
in this work two phenomenological dipole cross section 
parametrizations that include a realistic impact parame- 
ter dependence. One is the IIM [15] dipole cross section 
which is a parametrization including the most important 
features of BK evolution. The detailed expression for 
the dipole cross section can be found in Ref. [15]; we use 
here the values of the parameters from the newer fit to 
HERA data including charm from Ref. [16]. The second 
parametrization used here is a factorized approximation 
of the IPsat model with an eikonalized DGLAP-evolved 
gluon distribution [17, 18]. 

In the IIM model the impact parameter dependence is 
explicitly factorized as 
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(b T ,r T ,x) = 2T p (b T )AA(r T ,x), 



(2) 



We take, following Ref. [19], a Gaussian profile for the 
proton impact parameter profile function: T p (br) = 
exp (~b 2 /2B p ) with B p = 5.59 GeV~ 2 . 

In the IPsat model the impact parameter dependence 
is included in the saturation scale as 



da 



dip (b T , r T , x) = 2 [1 - exp (-r 2 F(x, r)T p (b T ))] 



d 2 b 

denoting r = |ry|. Here T p (b>r) is the same impact pa- 
rameter profile function as above, but the fitted value 
for the proton shape is B p = 4.0 GeV 2 (see Ref. [20] for 
a discussion about the different numerical value) and F 
is proportional to the DGLAP evolved gluon distribu- 
tion [21], 
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with C chosen as 4 and /i 2 , = 1.17 GeV resulting from 
the fit [18]. Following Ref. [20] we replace Eq. (3) by the 
factorized approximation 



d 2 b T 



(b T ,r T ,x)«2T p (b T ) [1 



exp(-r 2 ^(a;,r))] , 
(5) 

using the same F(x,r) defined in Eq. (4). This ap- 
proximation, denoted here as "flPsat", brings the IP- 
sat parametrization to the form Eq. (2) with Af(r, x) = 
[1 -cxp(-r 2 F(>,r))]. It was shown in Ref. [20] that 
also the flPsat parametrization also describes the HERA 
J/ty data accurately. 



III. DIFFRACTIVE CROSS SECTION IN 
ULTRAPERIPHERAL COLLISIONS 



allowed to break up. The event is still diffractive (there 
is a rapidity gap) as long as there is no exchange of color 
charge. 

The cross section for quasielastic (coher- 
ent +incoherent) vector meson production in nuclear DIS 
is (see e.g. [18]) 
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where — Q 2 is the virtuality of the photon. The coherent 
cross section is obtained by averaging the amplitude be- 
fore squaring it, | (A) N \ 2 , and the incoherent one is given 



by the variance (|„4| 2 ) N 
where 



(A) N | 2 , (see Refs. [20, 22]) 



r A 

<0({b Ti })) N = J JJ [d 2 b Tl T A (b Tl )] 0({b Ti }) (7) 



is the average over the positions of the nucleons in the 
nucleus. Here Ta is the Woods-Saxon distribution with 
nuclear radius R A = (1.12A 1 / 3 - 0.86A" 1 / 3 ) fm and sur- 
face thickness d — 0.54 fm. 

The factor 1 + (3 2 accounts for the real part of the 
scattering amplitude and the factor R 2 corrects for the 
skewedness effect, i.e. that the gluons in the target are 
probed at slightly different xp [23] . For these corrections 
we follow the prescription of Ref. [24] , taking them as 



tan ■ 



ttA 



Ra 



A = 



2 2A + 3 T(A + 5/2) 
A T(A + 4) 
din .4 



with 
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We calculate, as in Ref. [20], the correction terms from 
the energy dependence of the nucleon scattering ampli- 
tudes and use the same values for the nucleus at the 
same Q 2 ,xp. The real part and skewedness corrections, 
especially R g , are a significant factor in the absolute nor- 
malization of the cross section and are necessary for an 
agreement with HERA data. 

The imaginary part of the scattering amplitude, A, is 
the Fourier-transform of the dipole-target cross section 
(Tdip from impact parameter by to momentum transfer 
At, contracted with the overlap between the vector me- 
son and virtual photon wave functions: 



A(x P ,Q 2 ,A T ) = J d 2 v T J ^ J d 2 b 2 



x {^)(r,Q 2 lZ )e 



2 \ -ib-r-Ax 
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^(b T ,r T ,x P ), (11) 



In this work we consider both coherent and incoherent 
diffractive vector meson production. In a coherent pro- 
cess the nucleus off which the photon scatters remains 
intact, whereas in incoherent diffraction the nucleus is 



where we have followed the normalization convention of 
Ref. [18]. For the virtual photon-vector meson wavefunc- 
tion overlap we use the "boosted Gaussian" and "gaus- 
LC" parametrizations from Ref. [18]. 
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Assuming a large and smooth nucleus the averaged 
amplitude required to compute coherent J/ VP production 
reads [17] 

(A(x r ,Q 2 ,A T )) N = J ^d 2 r T d 2 b Te - b -^ 

x [*Jr*](r, Q 2 , z) 2 [1 - exp {~2irB p AT A (b)N(r, x P )}] . 

(12) 

At large —t = Ay 2 the cross section is almost purely 
incoherent. Thus the incoherent cross section can at large 
|i| be computed as the total quasielastic cross section, by 
first squaring and then averaging the amplitude. The 
result is derived e.g. in Ref. [20] and reads 



Pb+Pb -^J/f+Pb+Pb, /%n =2.76 TeV 
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(13) 



Following Ref. [25] we factorize the diffractive vector 
meson production cross section in nucleus-nucleus (or 
proton-nucleus) collisions to the product of the equiv- 
alent photon flux generated by one of the nuclei and the 
photon-nucleus cross section: 



AA-tJ/VA _ 



M- (14) 



Here a i A ^ J /^ A i s the diffractive photon-nucleus cross 
section, u = (My/2)e y is the energy of the photon in 
the collider frame and My and y are the vector meson 
mass and rapidity. The explicit expression for the photon 
flux n(uj) (integrated over the impact parameter of the 
AA-collision hj\ A > 2R A ) can be found in Ref. [25]. In 
nucleus-nucleus collisions both nuclei can act as a source 
of photons that scatter off the other nucleus: 



( ^ a A 1 A 2 ^J/^A 

dy 



(yW A Hy) + n A H-yW A2 (-y) 



T lA 2 



(15) 

In proton-nucleus collisions the photon flux generated 
by a nucleus is computed requiring that the impact pa- 
rameter is larger than R A . The proton can also act as a 
photon source, and the photon flux generated by a pro- 
ton is computed as in Ref. [25]. As the photon flux is 
proportional to the charge squared, the process where 
the photon is emitted from the nucleus dominates. 

The kinematics of diffractive vector meson production 
is such that the gluon xp probed by the real photon is 
At forward and backward rapidi- 



xp = M\ 

ties we have two different contributions: either a small- 
x photon scatters off a large-x gluon or vice versa. At 




FIG. 1: The coherent diffractive J/$ photoproduction (Q = 
GeV 2 ) cross section in lead-lead collisions at snn = 2.76 
TeV computed using flPsat and IIM parametrizations and 
Boosted Gaussian (thin blue lines) and Gaus-LC (thick black 
lines) wavefunctions compared with the ALICE data [10, 26]. 
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FIG. 2: The incoherent diffractive J/Sb photoproduction 
cross section in lead-lead collisions at ^/snn = 2.76 TeV com- 
puted using flPsat and IIM parametrizations and Boosted 
Gaussian (thin blue lines) and Gaus-LC (thick black lines) 
wavefunctions. 



midrapidity we only probe small- x structure of the nu- 
cleus. Our results should be most reliable in that region. 
At the LHC ^snn = 2.76 TeV, and for J/W production 
x P w 0.001 at y = 0. 
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FIG. 3: The coherent (thick lines) and incoherent (thin lines) 
diffractivc J/\& photoproduction cross section in lead-lead col- 
lision at ^/snn = 2.76 TeV as a function of momentum trans- 
fer t at midrapidity y = using the Gaus-LC wave function. 
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FIG. 4: The diffractive J photoproduction cross section 
in proton-lead collisions at ^/snn = 5.02 TeV computed using 
flPsat and IIM parametrizations and Boosted Gaussian (thin 
blue lines) and Gaus-LC (thick black lines) wavefunctions. 
The proton is moving in the negative y direction. 



IV. RESULTS AND DISCUSSION 

The ALICE collaboration has measured the coherent 
J/ty photoproduction cross section da/ dy at a relatively 
large rapidity \y\ ~ 3 [10]. Recently preliminary results 
at midrapidity y — were also published [26] . The com- 
parison between our results and the ALICE data is shown 
in Fig. 1. The PHENIX collaboration has also measured 
the ultraperipheral J/^F cross section in gold-gold colli- 
sions at ^/snn = 200 GeV and y = 0. For these kinemat- 



FIG. 5: Incoherent diffractive photoproduction cross 
section in lead-lead collision normalized by coherent J/vp pro- 
duction in pA (nuclear transparency ratio) at ^/snn = 2.76 
TeV computed using flPsat and IIM parametrizations and 
Boosted Gaussian (thin blue lines) and Gaus-LC (thick black 
lines) wavefunctions. The proton is moving in negative y di- 
rection. 



ics the flPsat dipole cross section with the gaus-LC wave- 
function gives the result 109/xb, compared to 76 ± 34/ib 
measured by PHENIX [27]. The ALICE data seems to 
favor the flPsat over the IIM model. This is perhaps not 
surprising. The most important difference between the 
two dipole models for this purpose is the different impact 
parameter dependence. The IIM value B p = 5.59 GeV -2 
comes from a fit to inclusive data and is close to the mea- 
sured value for inclusive diffraction. The HERA data 
for diffractive J/^> production [28, 29] has, however, a 
smaller B p , which is reflected in the IPsat parametriza- 
tion. Thus we consider results obtained using the flPsat 
model more reliable. 

Our results slightly overshoot the data, but the rapid- 
ity dependence comes out correctly. For example the flP- 
sat model with Gaus-LC wavefunction is above all data 
points by a factor ~ 1.4. We consider the agreement 
relatively good given the simplicity of the parametriza- 
tions and the fact that no nuclear data was used to con- 
strain the models. We emphasize that the parametriza- 
tions used are exactly the same as in Ref. [20], which 
predates the ALICE data. Recall that both dipole mod- 
els and wavefunctions used here give good descriptions 
of the HERA data. We suspect that the main reason for 
the larger normalization lies in the skewedness correction. 
The correction is larger here than for HERA kinematics 
because of the larger x probed, making it less reliable. 

The difference between the two wavefunctions is largest 
at Q 2 = which is the case here. The ALICE data 
seems to favor the Gaus-LC wavefunction, and thus we 
consider the flPsat dipole model and the Gaus-LC wave- 
function to be the most reliable combination. At \y\ > 2, 
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^/snn = 2.76 TeV, we are probing gluons with xp > 0.01, 
and in that region our parametrizations for the dipole 
amplitude are not valid any more. In addition, the 
real part and skewedness corrections in total become of 
the order 2, making them less reliable. Nevertheless, 
all wavefunctions and dipole models consistently give 
dcr/dy\ y=0 / da / dy\ y=2 = 1.41 ... 1.46. Thus the pre- 
diction for the rapidity dependence is much more robust 
than for the absolute normalization. 

We then present our predictions for the incoherent 
diffractive vector meson cross section in Fig. 2. Again 
the different models give a quite different overall normal- 
ization, but a very similar rapidity dependence. Notice 
that now the absolute normalization is larger in the flP- 
sat model. This is due to the different different impact 
parameter profiles, we refer the reader to Ref. [20] for 
a more detailed discussion. Again most of the difference 
cancels in the ratio da/ dy\ y =o / da/ dy\ y= z which is now 
1.35 . . . 1.43, so the energy dependence is very similar in 
coherent and incoherent scattering. 

Our result for the incoherent vector meson production 
from Ref. [20] is not valid at small \t\. However, we ex- 
pect to get a realistic estimate for the total incoherent 
cross section by integrating the differential cross section 
starting from the value of \t\ where incoherent and co- 
herent cross sections are equal. The error made is small, 
parametrically a factor ~ g-ltmmlBp w j^h \tmxa\ ~ V-^vl 
and numerically < 10%. 

In Fig. 3 we present predictions for the t distribution 
of diffractive J / \& photoproduction at midrapidity where 
xp 0.001. Note that the i-slope of the incoherent cross 
section directly measures the spatial distribution of glu- 
ons inside a nucleon, because the ^-dependence of the 
incoherent cross section is ~ exp(B p t) [20]. 

Our pA results arc shown in Figs. 4 and 5. In Fig. 4 
we show the rapidity dependence of the diffractive J/ 1 ^ 
cross section (the photon-nucleus scattering is required 
to be coherent). Now the difference between the mod- 
els is reduced as the dominant process is photon-proton 
scattering where the models are constrained by HERA 
data. Finally, in Fig. 5 we compute the incoherent J/^ 



photoproduction cross section in AA collisions divided 
by A times the diffractive J/^> production cross section 
in pA collisions at the same ^snn- Since ultraperiph- 
eral proton-nucleus collisions are mostly photon-proton 
collisions, this is a "nuclear transparency" ratio that 
measures the absorption of the dipole as it propagates 
through the nucleus, see e.g. [20, 30]. 



V. CONCLUSIONS 

We have computed coherent and incoherent diffractive 
J/$! photoproduction cross sections in ultraperipheral 
heavy ion collisions. The only inputs to our calculation 
come from fits to HERA data and standard nuclear geom- 
etry. Especially the rapidity dependence agrees relatively 
well with the published ALICE result, considering the 
rather large dependence on the details of the dipole model 
and the vector meson light cone wavefunction. The nor- 
malization of the data favors the flPsat parametrization, 
where the proton diffractive slope is constrained by the 
HERA diffractive J/\P data. We also find that the dif- 
ferent parametrizations, each known to fit HERA data 
well, yield significantly different normalizations for the 
cross section leaving the overall rapidity dependence very 
similar, with da/ dy\ v =o / da/ dy\ y =2 = 1-41 . . . 1.46 for 
the coherent and da j dy\ y= § / da/ dy\ y= 2 — 1.35 ... 1.43 
for the incoherent scattering. A similar conclusion was 
found in Ref. [20] for photon-nucleus scattering. We also 
present predictions for diffractive J /^ production in pA 
collisions. 
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